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Abstract 

The main purpose of this paper is to investigate the exact solutions 
of cylindrically symmetric spacetime in the context of f(R,T ) gravity 
PQ, where f(R,T ) is an arbitrary function of Ricci scalar R and trace 
of the energy momentum tensor T. We explore the exact solutions 
for two different classes of f(R,T) models. The first class f(R,T) = 
R + 2f(T) yields a solution which corresponds to an exterior metric of 
cosmic string while the second class f(R, T) = fi(R) + /RT) provides 
an additional solution representing a non-null electromagnetic field. 
The energy densities and corresponding functions for f(R,T ) models 
are evaluated in each case. 
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1 Introduction 

The most popular phenomenon in the modern day cosmology is the current 
expansion of universe. Observational and theoretical facts suggest that our 
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universe is in the phase of accelerated expansion [2]. The existence of dark 
matter and dark energy is another interesting topic of discussion [3j. Almost 
a century ago, Einstein gave the concept of dark energy by introducing a 
small positive cosmological constant in the held equations. But he rejected 
this idea later on. However, it is now believed that the cosmological con¬ 
stant may become a suitable candidate for dark energy. Modified theories 
of gravity seem attractive to explain late time acceleration of the universe. 
An interesting modified theory of gravity is the f(R) theory which involves 
a generic function of Ricci scalar in standard Einstein-Hilbert Lagrangian. 

In recent years, f(R ) gravity has been investigated by many authors in 
different contexts bhu]. Some interesting review articles [LS] can be helpful 
to understand the theory. Bamba et al. [US] explored curvature singularity 
appearing in the collapse process of a star in this theory. They established 
that curvature singularity could be avoided by adding R a term in the viable 
f(R) gravity models. Thermodynamics of the apparent horizon in the Pala¬ 
tini formalism of f(R) gravity has been discussed by Bamba and Geng [20J . 
Capozziello et al. [2Tj used Noether symmetries to find spherically symmetric 
solutions in f(R) gravity. Cylindrically symmetric vacuum and non-vacuum 
solutions have also been investigated in this theory [22] , Sharif and Shamir 
[ 23 ] explored plane symmetric solutions in metric f(R) gravity. The same 
authors [ 23 ] found the solutions of Bianchi types / and V cosmologies for 
vacuum and non-vacuum cases. Kucukakca and Camci [[22] discussed Pala¬ 
tini f(R) gravity using Noether gauge symmetry approach. For this purpose, 
they considered a flat Friedmann-Robertson-Walker (FRW) universe and it 
was concluded that the resulting form of f(R) model yielded a power law 
expansion for the scale factor of the universe. Conserved quantities in metric 
f(R) gravity via Noether symmetry approach have been calculated recently 

m- 

Another modified theory known as f(R, T) gravity has been developed by 
Harko et al. [T] . In fact, it is the generalization of f(R) theory of gravity and 
based upon the coupling of matter and geometry. In this theory, gravitational 
Lagrangian involves an arbitrary function of the scalar curvature R and the 
trace of the energy momentum tensor T. The equations of motions after the 
addition of an appropriate function f(T) indicate the presence of an extra 
force acting on test particles. The investigation of perihelion shift of Mercury 
using f(R,T) gravity provide an upper limit on the magnitude of the extra 
acceleration in the solar system which indicates the presence of dark energy 
[1]. Thus the study of f(R, T) gravity models may also provide better results 
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as compared to the predictions of standard theory of general relativity (GR). 
The action for f(R,T) theory of gravity is given by [1] 


s= f T ) + L m )d*x, (1) 

where g is the determinant of the metric tensor g lw and L m is the usual 
matter Lagrangian. It would be worthwhile to mention that if we replace 
f(R,T ) with f(R), we get the action for f(R) gravity and replacement of 
f(R, T ) with R leads to the action of GR. The energy momentum tensor T IJV 
is defined as [27] 

rp _ _2 8{y/^gL m ) 

V~9 1 j 

When we assume that the dependance of matter Lagrangian is merely on the 
metric tensor g IJV rather than its derivatives, we get 

= L m g lMJ - (3) 


Many authors have investigated this theory in recent years and a reasonable 
amount of work has been done so far. 

Adhav [28] explored the exact solutions of f(R, T) field equations for 
locally rotationally symmetric Bianchi type I spacetime. Bianchi Type V 
cosmology with cosmological constant has been studied in this theory by 
Ahmed and Pradhan [29]. Jamil et al. [30] reconstructed cosmological 
models in f(R,T) gravity and it was concluded that the dust fluid re¬ 
produced ACDM, phantom-non-phantom era and the phantom cosmology. 
Godel type universe was studied in the framework of f(R, T) gravity by San¬ 
tos pop. Sharif and Zubair [32] discussed the reconstruction and stability 
of f(R,T) gravity with Ricci and modified Ricci dark energy. The same 
authors [33] analyzed the laws of thermodynamics in this theory. However, 
it has been proved that the first law of black bole thermodynamics is vio¬ 
lated for f(R,T ) gravity [Mj. Houndjo [35] reconstructed f(R,T) gravity 
by taking f(R, T) = f\(R.) + f 2 (T ) where it was shown that f(R, T ) gravity 
allowed transition of matter from dominated phase to an acceleration phase. 
In a recent paper [36] , Harko and Lake investigated cylindrically symmet¬ 
ric interior string like solutions in f(R,L m ) theory of gravity. We explored 
Bianchi type / cosmology in f(R,T) gravity with some interesting results 
[37]. It was concluded that equation of state parameter w —» — last —> oo 
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which suggested an accelerated expansion of the universe. Thus it is hoped 
that f(R,T) gravity may explain the resent phase of cosmic acceleration 
of our universe. This theory can be used to explore many issues and may 
provide some satisfactory results. 

In this paper, we are focussed to find the exact solutions of cylindrically 
symmetric spacetime in the framework of f(R, T) gravity. The plan of paper 
is as follows: In section 2 , we give some basics of f(R,T) gravity. Section 3 
provides the exact solutions for cylindrically symmetric spacetime using two 
different classes of f(R,T) models. Summary and concluding remarks are 
given in the last section. 


2 Some Basics of f(R, T) Gravity 

The f(R, T) gravity field equations are obtained by varying the action S in 
Eq. (PP) with respect to the metric tensor g IIU 


fn(R, T)R^--f(R, T ) = KT^-f T (R, T)(T^+0^), 


( 4 ) 


where denotes the covariant derivative and 


df R {R,T) 


□ ee V"V M , f R (R, T) = fr(R, T) = 0^ = g afi —^ 


dfn{R,T) 


STa 

5g> w 


dR ’ J v ' dT 

Contraction of Eq.flJJ) yields 

fn(R, T)R + 3 Df R (R, T ) - 2 f{R, T) = kT - f T (R, T)(T + 0), (5) 


where 0 = 0^^. This is an important equation because it provides a rela¬ 
tionship between Ricci scalar R and the trace T of energy momentum tensor. 
Using matter Lagrangian L m , the standard matter energy-momentum tensor 
is derived as 

= 0 + v)u^u v - pg^u, ( 6 ) 

where u^ = yd?oo(l, 0, 0, 0) is the four-velocity in co-moving coordinates and 
p and p denote energy density and pressure of the fluid respectively. Perfect 
fluids problems involving energy density and pressure are not any easy task 
to deal with. Moreover, there does not exist any unique definition for matter 
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Lagrangian. Thus we can assume the matter Lagrangian as L m = —p which 
gives 

= ~P9nv ~ 2 T^, (7) 

and consequently the field equations (J4j) take the form 

fn(R, T)R^ u ~f(R 1 T)g^-(S7^v-g^^f^R, T) = KT^+f T (R, T^T^+pg^), 

( 8 ) 

It is mentioned here that these field equations depend on the physical nature 
of matter field. Many theoretical models corresponding to different matter 
contributions for f(R, T) gravity are possible. However, Harko et al. [T] gave 
three classes of these models 

r R+mT), 

f{R,T)=l fi(R) + h(T)i 
( h{R) + f2{R)MT). 

In this paper we are focussed to the first and second class. 


3 Exact Cylindrically Symmetric Solutions 

The line element of cylindrically symmetric spacetime is given by [38* [39] 

ds 2 = Adt 2 — dr 2 — B(dd 2 + a 2 dz 2 ), (9) 

where A and B are functions of radial coordinate r and a is an arbitrary 
constant. It may be pointed out here that Azadi et al. 122! have used 
another definition of cylindrical symmetry which may prove to be useful in 
some situations, but for our purpose we keep the above definition. Further, 
such type of line elements with cylindrical symmetry describe the spacetimes 
of a cosmic string. The corresponding Ricci scalar is 

1 T 2 A" A' 2 2 A'B' 4 B" B' 2 

~ 2 [~A~ ~ A 2 + AB + ~B ~ H 2 

where prime denotes derivative with respect to r. The main source of gravita¬ 
tional field is the energy-momentum tensor. For an ordinary star possessing 
cylindrical symmetry, p « p |40| . Therefore we can neglect pressure to 


( 10 ) 
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solve highly non-linear differential equations in this theory. Thus the energy- 
momentum tensor for dust is 


f 'jiw pU^U u , (ll) 

where p is the matter density and the four velocity vector u M satisfies the 
equation u M = 5 ° M . Now, we explore the solutions of the field equations for 
two classes of f(R,T ) models. 


3.1 f(R,T) = R + 2f(T) 

For the model f(R, T) = R + 2f(T), the field equations become 


R 


fllP 



+ 2f T (T)T fiU + 


f(T) + 2pf T (T) 


9ru- 


( 12 ) 


Here we find the most basic possible solution of this theory due to the com¬ 
plicated nature of held equations. However, in the next subsection we will 
investigate the solutions with more general case. For the sake of simplicity, 
we use natural system of units (G — c — 1) and f(T) = XT, where A is an 
arbitrary constant. In the case of dust with p — 0, the gravitational held 
equations take the form 

R/xv — + 2AT)g M! , — (8vr + 2A )T M „. (13) 


Thus for cylindrically symmetric spacetime, we obtain a set of differential 
equations for unknown A, B and p all depending on r. 


A'B' , 2 B" 
~AB + “B - 


B" 

B' 2 

= (87T + 3A )p, 

(14) 

~ ~B 

+ AB 2 

2 A'B' 

B' 2 

= -AXp, 

(15) 

AB 

+ ~B 2 

2 A" 

B' 2 


(16) 

+ A 


= —AXp. 


The conservation equation for energy momentum tensor is given by 


V lv . v = 0 , 


( 17 ) 
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which yields p = 0. Thus we obtain the metric coefficient A = constant 
as p 7 ^ 0. Without loss of generality, we take A = 1 and the two nonlinear 
differential equations (1151161) now reduce to 


2 B h 


B 


B' 2 

—4A p, 

(18) 

B' 2 

H 2 

—4A p. 

(19) 


Subtracting these equations, we get 


B' 2 B" _ 

~ ~B = °’ 


which yields a solution 


( 20 ) 

( 21 ) 


B(r) = 

where Ci and C 2 are integration constants. The metric takes the form 

ds 2 = dt 2 — dr 2 — cie C2r (d6 2 + a 2 dz 2 ), ( 22 ) 


This solution corresponds to an exterior metric of cosmic string [JT]. Energy 
density of universe p and trace of energy-momentum T in this case turn out 
to be 

o. 


3co 2 


p = 


4(8vr + 3A) 


= T, 


(23) 


while the Ricci scalar R becomes 


R = -C2 2 + 0. 


(24) 


The string type solutions in f(R, L m ) modified gravity also yields a constant 
Ricci scalar R = 0 [36]. In our case the Ricci scalar in non-zero constant, 
which is due to non-Kasner type nature of our solution. Here f(R,T) turns 
out to be 

f(R,T) = jfL ( 8,r + 2 \ (25) 

; 2 V 8 vr + 3A ; V ; 

and the mass of the string per unit length with radius r 0 can be calculated 

as (HI H 2 ] 


r2ir 


no 


m = 


d(f> / Ty^—gdr. 


(26) 
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Using Eqs. fllH I23|) . the mass per unit length of the string turns out to be 


37rc 1 c 2 2 ae C2r ° 
2(8vr + 3A) 


(27) 


We have to take A < — ^ to get positive energy density and mass of the 
string. It is mentioned here that when A = 0, this solution corresponds 
to constant curvature solution of cylindrically symmetric spacetime in f(R) 
gravity [43]. 


3.2 f(R,T) = fi(R) + h(T) 

Now we explore the solutions with more general class. Here we also take the 
dust case and the held equations for this model become 


= KT^+/ 2T (T)T^+i/ 2 (r)g„„. 

(28) 

Contracting the held equations, we obtain 

Rfin(R) - 2/i (R) + 3Df lR (R) = kT + Tf 2T (T) + 2 f 2 (T). (29) 


Using this, we can write 

, _ 3D fi R (R) + Rf inCg) -kT- TR t {T) - 2R(T) 

fi( ) 2 

Inserting this in Eq. (l28jh we get 

fin(R)R, u - ^ + f2 T {T))T^ 

RJir(R) - af 1R (R) -kT- T/ 2 t (T) - 4/ 2 (T) 


(30) 


(31) 


Since the metric ([9]) depends only on r, one can view Eq. (j3Tj) as the set of 
differential equations for /^(r), A and B. It follows from Eq. (j3TT) 

that the combination 


i _ fiR(R)Rfj,tJ, 
— 


\7,VJ lR (R)-(K + f 2T (T))T, u 
9rr 


( 32 ) 







is independent of the index /i and hence — A u = 0 for all /j, and u. Thus 
A 0 — A 1 = 0 gives 


A'B' A" B" B’ 2 A'fajR) f&(R) (« + / 2 r (r))p 

2 AB 4A B 2 B 2 2 AF F fm(R) 

Similarly A 0 — A 2 = 0 yields 


A'B' B" A’ 2 A" /^(R) A’ B' ( R + h T (T)),, 

4AB 2 B 4A 2 2A 2f ln (R) ' A I! Jir(R) ' 

Here we also assume the metric coefhcient A = constant due to the con¬ 
servation of energy-momentum tensor. Thus the two nonlinear differential 
equations (1331134|) now reduce to 


B" 5° /ifi(fl) (« + A r{T))p 

B 2 B 2 f 1R (R) 


B" B'h («+ JWCOjp 

2B 2B/ la (fl) 

Subtracting these equations, we obtain 

a' 2 B" 2/i«(fi) 

a 2 a Bf 1R (R) f IR (R) 


(36) 


(37) 


It has been established that dust matter-dark energy combined phases can 
be obtained by the exact solutions derived from a power law f(R) model |45j . 
Therefore, we follow the approach of Nojiri and Odintsov [46] and use the 
assumption fi R (R) oc foR m , where f 0 and m are arbitrary real constants. 
Thus Eq. (1371) takes the form 


2m(m 


lV R 2 

D w + m 



B'R'\ B" B’ 2 
BR ) + IF ~~ ~ER 


0 . 


(38) 


3.3 Exponential Solution 

For the constraint 

m(2m + 1) = 0, (39) 

Eq. (1381) admits the same exponential solution as given by Eq. (l2TT) . Two cases 
arise from the constraint. 
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Case I 


When m — 0, we obtain fi R (R ) = /o- This corresponds to fi(R) = / 0 -R + C3, 
where C 3 is an integration constant. For C 3 = 0, /o = 1 and / 2 (F) = 0, 
the solution correspond to GR. Here Ricci scalar is same as given in Eq. ([23j) 
while matter density and f(R, T) turn out to be 


~ c 2 (/qF + c 3 ) 

2(K + / 2T (r))’ 


f(R,T) = f 0 R + f 2 (T) + c 3 . 


(40) 


Many explicit expressions for energy density are possible with different choices 
of / 2 (F)- For example when / 2 (F) = XT 2 , it follows that 


K 


P~ + ^P + k = 0 , 


(41) 


where k = c P 2 'P , ^ +2c P . Two real roots of the quadratic equation (HID turn 
out to be 

-K ± vV - 16fcA 2 2 2 n 

p = -, k 2 - 16k\ 2 > 0. (42) 

4A 

Thus all the physical quantities turn out to be constant in this case. 


Case II 

We get f\ r {R) = foR— when m = — 1 / 2 . Here f\(R) becomes 


fi(R) = 2f 0 VR + C4, (43) 

where C 4 is an integration constant. The matter density and f(R,T) become 

/(fl.r) = 2/„v / l+/ 2 (T) + c 4 , (44) 

2(« + /2t(F)) 


P = 


This case also yields all the quantities constant. 


3.4 Power Law Solutions 

Here we assume the solution of Eq. ([38j) in power law form, i.e. B{r) = 
(c 5 r + c 6 ) n , where c 5 , c 6 are arbitrary constants and n is an integer. We get 
a solution for n — 2 with the constraint equations 

c 6 = 0, (45) 
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( 46 ) 


4m 2 + 4 m — 1 = 0. 

The solution metric takes the form 

ds 2 = dt 2 — dr 2 — c 2, r 2 (d9 2 + a 2 dz 2 ), (47) 

This solution matches to a spacetime which represents a non-null electro¬ 
magnetic field f39j. The Ricci scalar becomes R—\ which is non-constant. 
It would be worthwhile to mention here that exactly similar Ricci scalar has 
been obtained for Kasner type cylindrically symmetric solutions in f(R, L m ) 
gravity |36j. Here we also have two cases corresponding to two different roots 
of Eq. (|46|) . 


Case I 

The first root of Eq. tHRj) turns out to be m = ~ 1 t,^ 2 . For this, we have 

f 1R (R) = f„R (48) 

It has been shown that the terms with positive powers of the curvature sup¬ 
port the inflationary epoch m ■ After integrating Eq. (l48li and substituting 
the value of R , we obtain 


fi(R) — fo(~z) 2 + C7, 


(49) 


where /o = 7^75 and C 7 is an integration constant. Matter density turns out 


to be 


-2T(2-y/2)/ 0 


P = 


fo < 0 . 


(50) 


(^/^COr^’ 

Here we consider f 0 < 0 to get positive energy density. For special case when 
fiiT) = AT, the expression for energy density turns out to be 


P = 


-2=^(2-V2)/o 


(k + A )r l +^ 

while f(R,T) takes the form 

-2^[Ak(V2 - 1) + A(5\/2 - 6 )]/o 


(51) 


f(R,T) = 


(k + A)r 1+ ^ 


(52) 
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Case (i) Case (ii) 



Figure 1: Behavior of energy density versus radial coordinate for / 2 (F) = XT 
with A = 1, /o — —1 and hi = 8tt. 


The graphical behavior of energy density is shown in Fig(l). Also, when 
/ 2 (F) = XT 2 , we obtain 


2Ap T up T 


2 * (2- y/2)f 0 

rpl+y/2 


0 . 


(53) 


Fig. (2) shows the graphical behavior of p vs radial coordinate r for two roots 
of this equation. 


Case II 

For the second root m = ~ 1 ~ v ^ , we get fi R (R ) = foR 1 2 Vr . This case yields 

fi{R) = /o(^) 2 + c 8, (54) 

where /o = 1 j^ and Cs is an integration constant, ft is worthwhile to 
mention hare that negative power of curvature serves as an effective dark 
energy supporting the current cosmic acceleration 03- Matter density for 
this case takes the form 


P 


—2 2 (2 + y/2)f 0 

(«+ ' 


(55) 
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Case (i) 


Case (ii) 



Figure 2: Behavior of energy density versus radial coordinate for / 2 (F) = XT 2 
with A = 1, fo — —1 and k = 8n. 


Here we also assume / 0 to be negative in both cases to get physically accept¬ 
able results. The expressions for f(R,T) in special case / 2 (F) = XT takes 
the form 


—2 2 (2 + \/2)/o 

(k + A )r 1-v/ 2 ’ 


(56) 


while f(R,T) takes the form 


f(R,T) 


-2=^ [Ak(V2 + 1) + A(5\/2 + 6)]/ 0 

(k + A)r 1_v/ 2 


Also, when / 2 (F) = AF 2 , it follows that 


(57) 


2\p 2 + Acp -h 


2 2 (2 + a/2)/o 

r l—v/2 


0 . 


(58) 


The graphical behavior for two roots of this equation is shown in Fig. (2). 


4 Concluding Remarks 

This paper is devoted to explore the exact static cylindrically symmetric 
solutions in f(R,T ) gravity. To our knowledge, this is the first attempt to 
investigate cylindrically symmetric solutions in f(R, T) gravity. In this work, 
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we consider two classes of f(R,T) models. Moreover, we assume the dust 
case to find the solutions. First we take f(R,T) — R + 2 /(F). This case 
yields a solution which corresponds to exterior metric of a cosmic string. The 
Ricci scalar R , function of Ricci scalar f(R,T) and matter density p are all 
constant for this solution. The results are similar to the string type solutions 
in f(R,L m ) modified gravity yielding a constant Ricci scalar R = 0 [ 36 ] . 
However, in our case the Ricci scalar in non-zero constant, which is due to 
non-Kasner type nature of our solution. 

The second class with f(R, T) = fi(R) + /2(F) is the more general choice 
to explore the solutions. We assume fi R (R ) oc foR m , where /o and m are 
arbitrary real constants. The corresponding held equations are solved using 
exponential and power law forms of metric coefficient. Exponential solution 
is exactly similar as in the case of first class of f(R, T) model yielding all 
the quantities constant. However, the power law solution is similar to a 
spacetime representing a non-null electromagnetic held already available in 
the literature. So the physical relevance of these solutions is obvious. The 
function of Ricci scalar f \ (/?,) contains positive power of curvature in the 
hrst case while the second case corresponds to negative power of curvature. 
It is worth mentioning here that the terms with positive power of curvature 
support the inhationary epoch while the term with negative power of cur¬ 
vature serve as an effective dark energy which supports the current cosmic 
acceleration m- Moreover, we have discussed two choices for /2(F) with 
two cases each. When /2(F) = AF, the matter density in hrst case goes to 
zero as r approaches infinity. However, matter density increases with the 
increase in radial coordinate in the second case. For the case /2(F) = AF 2 , 
we get two graphs corresponding to two different roots for energy density. 
We discard the case where energy density is negative while the behavior of 
energy density for second root is similar as in the case /2(F) = AF. It would 
be worthwhile to mention here that when /2(F) = 0 , this class corresponds 
to f(R) gravity model and the results agree with [ 38 ]. Many other expres¬ 
sions for different energy densities can be reconstructed with choices of /2(F). 
Thus, it is hoped that such cylindrical solutions in the context of modified 
f(R, T) gravity may provide some interesting features of general relativistic 
strings and other topological defects that may have formed as a result of a 
phase transition in the early universe. 
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